Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 2 February 2008 (MN WF$t style file v2.2) 



Non-LTE dust nucleation in sub-saturated vapors 
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ABSTRACT 

We use the kinetic theory of nucleation to explore the properties of dust nucleation in 
sub-saturated vapors. Due to radiation losses, the sub-critical clusters have a smaller 
temperature compared to their vapor. This alters the dynamical balance between 
attachment and detachment of monomers, allowing for stable nucleation of grains in 
vapors that are sub-saturated for their temperature. We find this effect particularly 
important at low densities and in the absence of a strong background radiation field. 
We find new conditions for stable nucleation in the n—T phase diagram. The nucleation 
in the non-LTE regions is likely to be at much slower rate than in the super-saturated 
vapors. We evaluate the nucleation rate, warning the reader that it does depend on 
poorly substantiated properties of the macro-molecules assumed in the computation. 
On the other hand, the conditions for nucleation depend only on the properties of the 
large stable grains and are more robust. We finally point out that this mechanism may 
be relevant in the early universe as an initial dust pollution mechanism, since once 
the interstellar medium is polluted with dust, mantle growth is likely to be dominant 
over non-LTE nucleation in the diffuse medium. 

Key words: dust, extinction 



1 INTRODUCTION 



Dust particles are one of the fundamental components of 
. ,_| . the interstellar medium (ISM) and an ever-present worry 
for observers due to their opacity at optical and UV wave- 
5-1 1 lengths (Cardelli, Clayton & Mathis 1989). The ISM of the 
5^ , Milky Way is polluted by a mixture of grains made of a va- 
riety of materials, likely dominated by carbonaceous grains, 
silicates, and small PAHs particles (Mathis, Rumpl & Nord- 
sieck 1977; Weingartner & Draine 2001). The dust proper- 
ties are supposed to be the result of dust formation in the 
outflows of evolved stars (e.g. Salpeterl977; Stein & Soifer 
1983; Mathis 1990; Whittet 1992; Draine 2003, and refer- 
ences therein) and subsequent evolution, and eventual dis- 
solution, in the ISM, mainly as the effect of shock waves that 
destroy the grains through sputtering (Draine 1989; McKee 
1989; Edmunds 2001). Alternative dust production sites are 
supernova explosions (Kozasa, Hasegawa & Nomoto 1989, 
1991; Todini & Ferrara 2001; Nozawa et al. 2003; Schneider, 
Ferrara & Salvaterra 2004; Bianchi & Schneider 2007) and 
quasar outflows (Elvis, Marengo & Karovska, 2002). 

The theory of dust nucleation in astrophysics is heavily 
influenced by the theory of the nucleation of phase tran- 
sitions in super-saturated vapors (Becker & Doring 1935; 
Feder et al. 1966; Abraham 1974). The theory had mild suc- 
cess in reproducing nucleation rates, but is still controversial 
in many aspects, especially because it extrapolates the prop- 
erties of macroscopic bodies to clusters of few molecules and 



because it extends the thermodynamic approach to systems 
with a handful of particles. To add to these problems, astro- 
physical dust nucleation requires chemical reactions, since 
grains of materials that do not have a vapor state do nucleate 
(think, for example, to the nucleation of olivines from silicon 
oxides and metals; Draine 1979; Gail & Sedlmayr 1986). An 
alternative approach is the so-called kinetic theory, which 
describes nucleation as the result of attachment and detach- 
ment of monomers from a seed cluster of n particles (atoms, 
molecules or radicals; Nowakowski & Ruckenstein 1991ab). 

Both the thermodynamic and the kinetic nucleation 
theory have been developed in conditions of true equilib- 
rium, i.e., when the two phases have the same temperature. 
In the astrophysical scenario, however, the temperature of 
the dust grains can be sensibly lower than the temperature 
of the gas in which they are embedded due to efficient ra- 
diation cooling (e.g., Draine 1981). This would seem to be 
irrelevant to nucleation theory, since a vapor needs to be 
already nucleated in order to have grains that can be colder 
than the gas phase. Even a sub-saturated vapor, however, 
has a large number of unstable clusters that form by random 
association of monomers (and rapidly evaporate). In this pa- 
per we study the effect of cooling of these proto-clusters in a 
sub-saturated vapor, and the effect this has on the balance 
between attachment and detachment of monomers. Using 
the kinetic theory of nucleation, we find that even largely 
sub-saturated vapors can nucleate, provided they are not 
immersed in a strong radiation field. We compute, albeit 
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under some controversial assumptions, the non-LTE nucle- 
ation rate. We show that, even though it is not as large as 
in super-saturated vapors, it can produce dust grains at a 
rate that can reproduce the average dust grain density in 
the Milky Way over a timescale of several million years. In 
addition we show that non-LTE effects can increase the nu- 
cleation rate in super-saturated vapors. Non-LTE nucleation 
could therefore provide a slow channel for dust formation, 
in which dust is built over a relative long time in a slowly 
evolving region. Such an example could be the outflow from 
AGN nuclei (Elvis et al. 2002). Such evolution is different 
from the one envisaged in the classical dust factories - AGB 
star atmospheres and supernovae - where dust nucleation 
is rapid but short lived since the favorable conditions are 
rapidly lost. 

This paper is organized as follows: in § 2 we briefly 
review the classical kinetic theory of nucleation; in § 3 we 
compute the dust grain temperature and in § 4 we compute 
the new conditions for nucleation. In § 5 we consider the 
nucleation rate and discuss our results in S 6. 



2 CLASSICAL NUCLEATION 

In this section we review some basic concepts of the theory of 
nucleation at densities and temperatures typical of labora- 
tory experiments (in conditions of "true equilibrium"). Let 
us first consider a flat surface separating the vapor phase 
of a certain material with its condensed phase (either liquid 
or solid) . The equilibrium condition between the two phases 
can be described in either thermodynamic or kinetic terms. 

In thermodynamic terms, equilibrium implies that the 
temperature of the two phases are equal and that the chem- 
ical potential of the molecules in the two phases are equal 
(e.g., Vehkamaki 2006). In this paper we adopt the kinetic 
approach, since we will investigate nucleation out of thermal 
equilibrium. In the kinetic representation the equilibrium 
is dynamical and is obtained by equating the rate of va- 
por phase molecules that become attached to the condensed 
phase to the rate of molecules that return from the con- 
densed phase to the vapor. The first rate is relatively easy 
to quantify, as long as it can be assumed that the vapor 
molecule velocity distribution is MaxwelliarQ: 
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dn in , / kT 
= ksnxy 



dt dA 



2ivmx 



(1) 



where nx is the density of the molecules (or atoms) of the 
compound of interest in the vapor phase, k is the Boltzmann 
constant, T is the temperature and rnx the molecular mass. 
The parameter k s < 1 is the sticking coefficient and repre- 
sents the probability that the incoming monomer remains 
attached to the condensed phase rather than bounce and 
return to the vapor. There is no reason for k s to be con- 
stant and it does indeed depend on the temperature (see, 
e.g., Batista et al. 2005). For the sake of simplicity, we will 
consider here k s to be a constant. 

The rate of detachment of monomers from the con- 
densed phase is much harder - if at all possible - to compute 



We concentrate here on homogeneous nucleation in vapors with 
only inert species and the compound of interest. 
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Figure 1. The attachment and detachment of monomers from a 
water droplet at T = 293 K and saturation S = 4. The solid curve 
shows the attachment rate, while the dashed curve shows the 
detachment rate as a function of the droplet radius. At small radii, 
the detachment dominates and the droplet tends to evaporate. 
At large radii, the attachment dominates and the drop tends to 
grow in size. There exist a radius where the rates are equal and 
the droplet is stable. This is called the critical radius r cr , and 
the droplet, called critical cluster, contains a critical number of 
monomers. 



a-priori. It can be obtained by comparing the experimentally 
determined saturation density nx,aq{T) with the theoretical 
attachment rate from Eq.[T](see, e.g., Kashchiev 2000, here- 
inafter K00, and references therein): 



dn 



dt dA 



= k B n x ,cq{T) 



kT 
2Tvmx 



(2) 



Equation [5] which is formally identical to Eq. [T] be- 
comes interesting if we assume that the detachment rate 
from the condensed phase does not depend on the proper- 
ties of the vapor phase (pressure, temperature, saturation), 
This seems to be a good assumption, since the properties 
and dynamics of a condensed material should be rathher 
independent from the low-density gas surrounding it. How- 
ever, the detachment rate may depend slightly on the vapor 
conditions (pressure, temperature, density) through their ef- 
fect on the properties of the condensed material (density, 
conductivity, surface tension). We neglect here these depen- 
dencies since there is no established theory to model them. 
With this assumption, Eq. [5] give us the detachment rate as 
a function of the properties of the condensed phase only, as 
long as the saturation density as a function of temperature 
is known either theoretically or experimentally. 

So far we discussed the equilibrium between a flat sur- 
face and a vapor. Phase transitions cannot take place in 
a single bulk event, but must be realized through micro- 
scopical nucleation of the new phase, due to energetic con- 
straints (K00). Equation [2] is modified if the surface that 
separate the vapor from the condensed phase is not flat. Let 
us consider a spherical droplet of radius r. Both thermo- 
dynamic (Vehkamaki 2006) and kinetic consideration (K00) 
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show that the saturation density for a droplet is larger than 
for a flat surface. In the thermodynamic case, this is due to 
the surface energy component, while in the kinetic approach, 
the effect can be explained through the fact that the binding 
energy of a surface monomer is smaller than for a bulk one, 
and so it is easier to eject a monomer for a small droplet 
than for a flat surface (at the same temperature). As long 
as the number of molecules in the droplet is large enough 
to allow for the definition of a surface, it can be shown that 
(K00): 



(3) 
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nx, eq (r,T) = nx, eq (T)e^ 



where a is the surface tension of the condensed phase and 
Do the volume occupied by one molecule in the condensed 
phase. Since the rate of impacts of vapor molecules on the 
surface does not depend on the curvature, and assuming that 
the sticking coefficient does not either, we can easily find the 
detachment rate as a function of temperature as: 

dn out (r) 

-nnnr = k ° nx ^ {T) 



kT 
2irmx 



2<rrp 
o rkT 



(4) 



where we used the so-called capillary approximation, i.e., 
we assume that a does not depend on r. Even though 
not explicit, the ejection rate of Eq. [4] depends strongly 
on the grain temperature through the equilibrium density 
nx,eq(T). The equilibrium density depends exponentially on 
the temperature, rougly as nx,cq(T) <x e~ A/lkT , where A is 
a positive, material dependent, constant. 

Figure [T] compares the attachment and detachment 
rates of monomers for water droplets at a temperature 



293 K and at saturation S 



nx/nx,, 



4. We 



adopted a = 78 dyne/cm, vo = m,H 2 o/p = 3 x 10 cm , 
k s = 1, and we obtained the equilibrium pressure from the 
CHERIC (Chemical Engineering Research Information Cen- 
ter) web sitfl At small radii, the detachment rate is larger 
than the attachment rate and the droplet will tend to evap- 
orate, while at large radii the droplet tends to grow. At the 
critical radius r cr the two rates are equal and the droplet 
is in equilibrium. The droplet with this radius is also called 
critical cluster. It is easy to see that the problem of nucle- 
ation reduces to the difficulty of creating droplets big enough 
to be super-critical, since once the critical radius is passed, 
the droplet will grow into the condensed phase. 

The fact that a droplet (or nucleus) smaller than r cr 
tends to evaporate does not mean that it is impossible to nu- 
cleate in the absence of pre-existing super-critical droplets. 
It exists, in fact, the possibility that by statistical fluctua- 
tions a nucleus grows big enough to reach the critical radius 
and subsequently evolve in a stable droplet. For a super- 
saturated vapor, it can be shown that the stationary distri- 
bution of nuclei is given by (K00 and references therein): 



dn x nr= 2 at 



\3=2 



where n is the number of molecules in the droplet, 
2 dnin 



/„ = 4-7T r" 



dt dA 



(6) 
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Figure 2. Nucleation rate for water as a function of saturation for 
three different temperatures: T = 260 K (dot-dashed line); T = 
293 K (solid line), and T = 310 K (dashed line). Experimental 
results from Wolk & Stray (2001) for T = 260 K are shown with 
solid dots. 

is the attachment rate for the droplet with n = (47r/3) r 3 /no 
molecules, and 

2 dn ou t ,_v 
9n =^ r dTdA (7) 

is the detachment rate for the same droplet. Equation [5] is 
quite complex. We refer to K00 for a more detailed discus- 
sion. We here simply notice that since the detachment and 
attachment rates enter only through their ratio, the station- 
ary distribution does not depend on the unknown sticking 
coefficient k s . 

Another fundamental result of the kinetic approach is 
the nucleation rate J, i.e., the rate at which new nuclei of the 
condensed phase are created in the vapor. It reads (K00): 



J — nx fi 



\j=2 



and it now depends linearly on the sticking coefficient k s 
through /i. In Eq. [5]and|8j the infinite limits of summation 
can be substituted to summation up to twice the number 
of molecules in the critical cluster (as long as only the dis- 
tribution up to the critical cluster is of interest). Finally, 
Eq. [5] and [8] are valid in the nucleation stage, i.e., when the 
grain growth and destruction proceeds through attachment 
or detachment of monomers. In the coalescence stage, when 
grain-grain collisions are relevant, different equations must 
be adopted. 

Figure [2] shows the nucleation rate as a function of sat- 
uration S for water droplets in vapor. The radius of the crit- 
ical cluster is computed inverting equating the attachment 
and detachment rates (Eq. [T]and|4]) and reads: 



2cjvq 



(9) 



http : //www. cheric . org/research/kdb/hcprop/cmpsrch.php 



kTlnS 

The fact that the nucleation rate increases with saturation 
is a direct consequence of the decrease of the critical radius. 
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A smaller critical radius requires less statistically unfavored 
random associations to reach the stable configuration and 
therefore increases the nucleation rate. In Eq. [8] the values 
of fi are directly proportional to saturation, while the gi do 
not depend on S. This difference brings about, at the math- 
ematical level, the dependence of the nucleation rate on sat- 
uration. Three temperatures are shown in Fig. [2] T = 293 K 
(solid), T = 260 K (dash-dot) and T = 310 K (dashed). The 
solid dots show the result of a water nucleation experiment 
(Wolk & Stray 2001) at T = 260 K. The comparison of the 
theory with the data shows that there is a discrepancy of 
several orders of magnitude between data and experiment. 
This could be due to k s -C 1. However, such an explanation 
is unlikely for water droplets (see Batista et al. 2005). It is 
rather believed that the discrepancy is due to the failure of 
the capillary approximation for droplets with a very small 
n. A larger surface tension for small droplets is required to 
make the theory consistent with observations. When com- 
menting the quantitative results of the nucleation theory, we 
shall keep in mind that reality may differ by few orders of 
magnitude from the theory. 



3 THERMAL BALANCE 

All the theoretical framework discussed in the above section 
was developed in an attempt to describe natural and labora- 
tory evidence of nucleation (e.g. rain). In the astrophysical 
environment, conditions can be very different from the lab- 
oratory and the theory needs some adjustments. First we 
consider that, in general, astrophysical vapors are immersed 
in an overwhelming quantity of inert H atoms (or, at least, 
inert in the condensation process). Second, at astrophysical 
densities the dust particles are at a lower temperature than 
the gas phase. In this section, we compute the effect of the 
non thermal equilibrium between gas and droplets (here- 
inafter dust grains). There are several components affecting 
the grain temperature, the most important ones are heating 
by collisions with gas phase particles and heating/cooling 
through the interaction with the radiation field. 

We here consider two effects: heating by collisions with 
gas particles and cooling due to the emission of radiation. We 
neglect radiation heating, i.e. we assume the radiation field 
is diluted and unimportant (radiation is a sink of energy). 
The cooling due to radiation is given by: 

dE 
~dt 



= 4ivr 2 asBA{r,Ts)Ts 



(10) 



where ctsb = 5.67 x 10 -5 erg cm -2 K -4 s -1 is the Stefan- 
Boltzmann constant, A < 1 a factor that takes into account 
that the grain does not emit as a black body at all frequen- 
cies and T s is the temperature of the grain. A(r, T) is given 
by (Laor & Draine 1993): 



A(r,T) = 



J Q* bB (\)Bx(T) d\ 
J B X {T) dX 



(11) 



where Q a bs(A) is the absorption efficiency and B\{T) is the 
black-body spectrum. 

To compute the heating of the grain due to collisions, 
we consider the astrophysical scenario, where most of the 
collision between the grain and the gas particles do not lead 
to attachment because of the overwhelming abundance of H 
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Figure 3. Temperature of water droplets in a vapor as a function 
of the vapor temperature T g . The temperature depends on the 
droplet radius since the radiation emitted is proportional to the 
droplet opacity, which scales linearly with the radius. From top 
to bottom, the lines show droplets of radius: 10 — 3 , 3 X 10~ 3 , 
10 -2 , 3 X 10 -2 , 10 _1 and > 1/im. These radii correspond to 
3,10,30,100,300, and 3000 inter-molecular distances, respectively. 



over the grain monomers. We also neglect vibrational and 
rotational energy of the gas particles (again, we suppose to 
be dominated by H). Finally, we neglect the sputtering of 
monomers from the cluster that can result from the collision 
with a high speed gas particle. Such effect is known to be 
relevant in shock destruction of dust particles but is usually 
neglected in the classical theory of nucleation (K00). A de- 
tailed treatment of the effect requires a detailed description 
of the energy levels of the clusters, a formidable task that is 
beyond the scope of this paper. We obtain: 

§ = A*r> j~?f P {v)(\» mH v>-\kT s )dv (12) 



AmrW —n g kT g 



TT/imH 



(13) 



where p(v) is the Maxwellian distribution of velocities, ma 
is the hydrogen mass, and [i ~ 1.2 the mean atomic weight 
of the gas. The subscript s always refer to the solid (grain) 
phase, while the suffix s always refers to the gas (vapor) 
phase. 

The balance between losses and gains of energy gives 
us an implicit equation for the temperature of the grain 
as a function of the gas pressure and temperature. Assum- 
ing a single equilibrium temperature for small grains is an 
oversimplification (Guhathakurta & Draine 1989). A broad 
range of temperatures would likely increase the nucleation 
rate since the cloder-than-average grains could grow into 
stable grains even for conditions where the average temper- 
ature would not allow for nucleation. On the other hand, 
hotter-than-average grains would be destroyed, providing a 
balancing effect. Neglecting the temperature distribution we 
obtain: 
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Figure 4. Same as Fig. [T]but for an sub-saturated vapor. Nucle- 
ation is still possible due to the cooling of the droplet. A critical 
radius can be found in this case, analogously to the classical nu- 
cleation scenario. 



A(r,T s )Tt = 



2kT g 
irfirriH 



1/2 



C"SB j 



/ 



(14) 



Figure [3] shows the magnitude of the effect for wa- 
ter droplet in a hydrogen dominated gas of total density 
n g = 10 10 cm" 3 . For the largest size we used the approx- 
imation A = 1, while for small water droplets we adopted 
the approximation: 



A{r,T) 



rr r ; 



r < r r 



(15) 



where r r is the opacity per unit length in the medium. For 
pure water we have A(r) ~ 1200 r in the temperature range 
250 < T < 600 K (pure water absorption coefficients are 
taken from the Oregon Medical Laser Center-OMLC web 
paglj. 



4 KINETIC EQUILIBRIUM 

As noted by Draine (1981), the departure from thermal equi- 
librium can have important effects on the kinetic balance of 
the clusters. Even though the effect of the grain radiative 
cooling is not exceedingly large (see Fig. especially for 
the small grains, the detachment rates change exponentially 
with temperature. This implies that even a small temper- 
ature change can make the difference between a nucleating 
mixture and a non- nucleating one. An additional complica- 
tion is caused by the fact that while the grain temperature 
depends on the total gas density n g , the attachment rate de- 
pends on the partial gas density of the compound of interest 
nx = Hx n g , where we have defined the number abundance 
of the compound Hx- 

Figure [4] shows the effect of the temperature change 
in the dynamic equilibrium for water droplets in a vapor 
with temperature T g = 293 K, n g = 5 x 10 8 cm -3 , and 



detachment rate would be ~ 10 10 , about ten orders of mag- 
nitude larger than the attachment rate, and droplets would 
evaporate. Due to the approximately exponential depen- 
dence of the experimental saturation density with tempera- 
ture, the change in temperature of the droplets suppresses 
the detachment rate. The situation appears now qualita- 
tively similar to the nucleating condition in Fig. [T] At small 
radii, the detachment rate is larger than the attachment 
rate, and the droplets evaporate. A critical radius r CT does 
exist, where the two rates are identical and the droplet is in 
equilibrium. The critical radius is now given by a modified 
version of Eq. [9] since the gas and grain temperatures are 
different: 



2avo 
kT 3 In Z 



where 



(16) 



(17) 



http : //omlc . ogi . edu/spectra/water/ 



In the same way as in the LTE case, droplets larger 
than the critical radius grow due to the larger attachment 
rate with respect to the detachment rate. 

Figure [4] was computed solving numerically Eq. [14] for 
the temperature with A(r, T) obtained from Eq. [11] with 
water data taken from the OMLC web page (see footnote [3}. 
Again, k a — 1 was assumed. Relaxing this assumption would 
change the normalization of the two curves so that r cr would 
remain constant. 



5 NUCLEATION 

The nucleation rate for the sub-saturated nucleating vapor 
can be computed using Eq. [5] given the attachment rate for 
the gas temperature and the detachment rate computed at 
the grain temperature. Figure [5] shows an extension of Fig. [2] 
to much lower saturations (and therefore densities). When 
saturation drops below unity, nucleation does not take place, 
until at very low densities and saturations non-LTE nucle- 
ation sets in. The nucleation rate is somewhat small, since 
the attachment rate is low at low saturation, but the detach- 
ment rate is suppressed by the decreased temperature of the 
grains and therefore nucleation is possible. It must be em- 
phasized that the theory becomes more and more inaccurate 
as we approach the limit of a small number of monomers in 
the critical cluster. In this case, many of the approximations 
made to derive the theory lose validity. First, the opacity per 
unit length looses meaning, and the cooling rate is not com- 
puted accurately. Second, the cluster does not have a well- 
defined surface any more and the capillary approximation 
for the surface tension breaks down. For this reason we did 
not include in Fig.[S]the very low saturation limit, since the 
number of monomers in the critical cluster becomes close to 
unity. Nucleation is still possible there, but computing the 
rate in the atomistic approximation is beyond the goal of 

4 The saturation is defined as S = nx /nx,eq(T g ) also in non- 
LTE conditions, but S = 1 loses in this case the meaning of 
threshold for the nucleation process. 
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Figure 5. Same as Fig. [2] but extended to very low saturation. 
The nucleation rate is non-negligible and nucleation is possible. 
The inset shows the number of monomers in the critical cluster as 
a function of saturation. Contrary to the classical case, the critical 
radius increase with the saturation, in this case. The curves on 
the right side of the plot are the results for LTE nucleation from 
Fig.fJ 



this paper and, for what matters, a very controversial issue. 
We show in the inset of Fig. [5] the number of monomers in 
the critical cluster as a function of saturation in the non- 
LTE nucleating region. Calculation have been stopped at 
the limit n C i = 2, but for n cr 5 the results should be taken 
with caution. 

Finally, we show in Fig. [6] the new conditions for water 
droplet nucleation in the n — T plane. As explained above, 
nucleation in the non-LTE area depends on the abundances 
Efr 2 o- In dark gray we show the region where nucleation is 
possible in the classical view. In light gray, we show the area 
where non-LTE nucleation is possible in the limit ~S.h 2 o — * 1- 
Solid lines show how the area is modified in the cases of 
Sh 2 o = 10 _1 , 10~ 3 , and 10 -5 . The three thin vertical lines 
in the lower left part of the diagram show the conditions 
under which nucleation rates have been computed in Fig. [5] 
It is important to note that even though nucleation is pos- 
sible at higher densities, the nucleation rate becomes non- 
negligible only for densities roughly equal or smaller than 
nn 2 o = 10 s cm -3 . The dependence of nucleation on density, 
for a given temperature, is quite complex. At very high den- 
sities, the conditions of temperature equilibrium are main- 
tained by the high rate of collisional heating and classical 
LTE nucleation takes place. As the density is decreased, two 
different things can happen. At high temperature (the re- 
gion in Fig. and [7] where a white area is present), the 
vapor becomes sub-saturated and nucleation is halted. De- 
creasing the density even further allows for the cooling of 
the sub-critical grains, and nucleation in non-LTE sets in. 
If the temperature is lower, non-LTE conditions apply even 
in saturated vapors. In that case, the vapor is always nu- 
cleating (see Fig. [SJ. The dashed line in Fig. [5] shows the 
region where the temperature of the droplets is equal to the 
freezing point of water (assuming it freezes at 0° C indepen- 
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Figure 6. Nucleation phase diagram for water droplets. The clas- 
sical vapor pressure line is shown as a thin solid line. The new 
equilibrium is shown as a thick solid line. Only the region inside 
the line does not produce any nucleation. The three thin vertical 
lines in the lower left corner show the values for which the nu- 
cleation rate is shown in Fig. \E\ The dashed line shows the ice 
formation boundary: only the droplets in the upper right corner 
are liquid (assuming that water freezes at 273.15 K independently 
of pressure). 

dent on pressure). Non-LTE nucleation results therefore in 
ice and not liquid droplets. Finally, the above treatment is 
valid only if there are no pre-existing grains (e.g. graphite 
or silicates) that can work as seeds for heterogeneous nu- 
cleation. In that case, the presence of already cold grains 
could result in the growth of ice mantels rather than in the 
nucleation of new droplets. 

5.1 Graphite 

Even though water has an astrophysical relevance for dust, 
especially in covering refractory grains with ice mantles, we 
now consider a more pregnant example: graphite grains con- 
densing from carbon atoms in the gas phase. For simplicity, 
we consider a gas made only of hydrogen with solar carbon 
abundance: E c = 3.3 x 10" 4 (Grevesse & Sauval 1998). We 
use the vapor pressure measurement from Brewer, Gilles & 
Jenkins (1948) and the surface tension ac = 34.6 dyne/cm 
(Morcos 1972). 

Figure [7] is analogous to Fig. [6] but shows the phase 
diagram of graphite. Note also that the y-axis reports the 
total density of the gas and not the partial density of carbon 
atoms. Analogously to the case of water, we see that graphite 
grains can form in non-LTE at temperatures much larger 
than in the classical scenario. 

To compute nucleation rates in non-LTE, we need to 
derive the factor A for graphite. We use the tables^ of Q 

abs 

from Draine & Lee (1984) and Laor & Draine (1993). In 
some cases, however, the grain sizes given in the tables are 

° http : //www . astro . princeton . edu/~draine/dust/dust . diel . html 
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Figure 7. Nucleation phase diagram for graphite grain in a hy- 
drogen gas with carbon number abundance E.q = 3.3 X 10 — 4 . 
As for the case of water analyzed above, there exist a big region 
at high temperatures and low densities where graphite grains can 
nucleate in non-LTE conditions with the gas. The point marks the 
conditions for which the nucleation rate J ss 7 X 10 4 cm -3 s _1 
has been computed in the text. 

not small enough (the smallest grain considered in the tables 
has r — 0.001/im while the smallest we consider has r = 
0.00013/im). In that case we extrapolate to smaller radii in 
the optically thin limit assuming Q a bs(A, r) = 1 — exp(— rr r ). 

In the case of graphite, nucleation in the non-LTE re- 
gion can be very efficient. For a gas density n 9 = 2.5 x 
10 13 cm -3 and temperature T g = 4500 K (see solid dot in 
Fig. [7}, the critical cluster has ~ 7 carbon atoms and the 
nucleation rate is J ~ 7 x 10 4 cm -3 s _1 . 

Another important aspect of Fig.[7]is that there is a vast 
temperature range where there is no non-nucleating region 
between the classical and non-LTE nucleating regions. This 
was true for water as well, but in a narrower region. For 
T g < 3250 K in the carbon-hydrogen mixture considered, 
the grain temperature is smaller than the gas temperature 
even in the classical nucleating region. This implies that 
the nucleation rate is different from the classical value and 
the cooling of the grains has to be considered even in the 
classical region. Figure [5] shows the nucleation rate for T = 
2500 K in the saturation range 0.1 < S < 10. Even though 
a mixture in true equilibrium does nucleate for 5* > 1, non- 
LTE effects increase the rate. 



6 DISCUSSION AND CONCLUSIONS 

We have considered the effect of grain (cluster, droplet) cool- 
ing in the nucleation of liquid and solid phases in vapors. We 
find that the effect can be dramatic on the nucleation rate 
and on the nucleation phase diagram, allowing for nucleation 
in large regions of the parameter space that are classically 
considered to be non-nucleating. 

As is in general true for nucleation, there are several lim- 
its and approximations that we should bear in mind when 



Figure 8. Nucleation rate in a carbon-hydrogen mixture with 
solar carbon abundance at T = 2570 K. The dashed line shows 
the classical result for "true equilibrium" nucleation, while the 
solid line shows the result of non-LTE nucleation. Even for S 1, 
where LTE nucleation takes place, allowing for the grain cool- 
ing increases the nucleation rate. The inset shows the number of 
monomers in the critical cluster for the two cases as a function 
fo the saturation. Not that computations in the super-saturated 
regime suffer the uncertainties of the few monomers limit for both 
LTE and non-LTE assumptions. 

considering the theory from the quantitative point of view. 
As exemplified by the comparison of the data with the the- 
ory in Fig. [2] several orders of magnitude can separate the 
nucleation rate prediction from the observations. The con- 
troversial points are: 

• Sticking coefficients — Most of the figures and compu- 
tations in this work assume k s — 1. This is not always true 
(Batista et al. 2005). In addition to its dependence on tem- 
perature, the sticking coefficient may depend on the size of 
the cluster. A big cluster could more easily absorb the ex- 
tra kinetic energy of the incoming monomer, compared to 
a small cluster (K00) , and therefore k s may be significantly 
smaller than unity for very small clusters. 

• Capillary approximation — The capillary approxima- 
tion, i.e., the assumption that the surface tension does not 
depend on the cluster size, is very controversial, and a 
change in the surface energy for very small clusters could 
result in big changes on the nucleation rate. For example, 
the discrepancy in Fig. [2] could be solved by assuming a 
larger surface tension for the very small water droplet. In 
addition, macromolecules do not even have a properly de- 
fined surface, and the whole concept does not apply. Finally, 
the exponential factor in Eq. [3] depends on the assumption 
that the clusters are spherical. A different ratio of the sur- 
face to the volume would modify this term. This is likely for 
small graphite clusters, since graphite tends to aggregate in 
a planar form. 

• Detachment rate for small clusters — In this paper, and 
in most nucleation theory, the detachment rate is computed 
by propagating to very small clusters the detachment rate of 
macroscopic bodies. It is very likely that, in the very small 
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limit, the detachment is governed by completely different 
processes. Let us analyze the two limits. For a macroscopic 
body, the number of monomers is so large that monomers 
with a statistically higher energy can detach since their en- 
ergy is larger than the binding energy. In the opposite limit 
of a dimer, the detachment has to be due to an external ac- 
tion: either a collision with a fast monomer or with a photon 
or with another cluster. In this limit destructive collisions 
have to be considered. 

• Cooling of macromolecules — This is a new problem 
that arises when the cooling of the grains is considered. We 
have assumed that the grains cool as modified black bodies 
down to the smallest sizes. When the number of monomers 
in the cluster becomes small, the cooling will not be through 
a continuum spectrum, but through lines and bands. A more 
refined treatment of cooling is necessary to compute accu- 
rate nucleation rates. 

• It should finally be kept in mind that we allowed for 
the complete cooling of the grains, neglecting the effects 
of a background radiation field in setting a lower limit to 
the temperature. In addition, we neglected the fact that the 
temperature of small grains is largely stochastic and that 
at high temperature some collisions between grains and gas 
particles can result in sputtering rather than accretion. 



Despite all these caveats, the main result of this pa- 
per holds: the region where a vapor spontaneously nucleate 
is not limited to the classic region, when nucleation takes 
place in thermal equilibrium and T g =T S . Allowing for the 
clusters to cool inhibits the detachment of monomers from 
the clusters and allow for nucleation at higher temperatures 
and lower densities than in the classical scenario. The nucle- 
ation rates tend to be orders of magnitude smaller than those 
derived in thermal equilibrium, but provide a non-negligible 
correction to the LTE rate for mildly super-saturated vapors 
with S ^ 1, especially at low temperature. 

In the astrophysical scenario small nucleation rates are 
not a big worry. The ISM of our Galaxy contains approxi- 
mately 1 per cent of its mass in dust grains (Mathis et al. 
1977). This corresponds to approximately one dust particle 
every cubic meter (assuming a power-law grain size distri- 
bution as in Mathis et al. 1977). However, in the present day 
Universe the ISM is already polluted with dust and in the 
presence of seed grains it is likely that the process of man- 
tle growth dominates over non-LTE nucleation in the diffuse 
medium. The process of non-LTE nucleation may therefore 
be important at high redshift, when the ISM is first pol- 
luted with metals by supernova explosions. It is unclear if 
supernovae do generate dust by themselves, and even more 
whether the generated dust can survive the sputtering in the 
forward-reverse shock systems (Bianchi & Schneider 2007; 
Nath, Laskar & Shull 2007). In the case that supernovas do 
mainly pollute the ISM with metals but with no or a negligi- 
ble quantity of dust, non-LTE nucleation could become the 
dominant process of dust nucleation in the early universe. 
Detailed estimates of nucleation in the various scenarios re- 
quire a more detailed understanding of the properties of the 
very small nuclei and are beyond the scope of this paper. 
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